Astronomy & Astrophysics manuscript no. 1694 


February 2, 2008 


(DOI: will be inserted by hand later) 





2 



The distance to the Pleiades 

Main sequence fitting in the near infrared 

S.M. Percival^, M. Salaris"*^, and M.A.T. Groenewegen^ 

Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, 
Birkenhead, CH41 ILD, UK 
e-mail: smp,ms@astro . livjm. ac.uk 

Instituut voor Sterrenkunde, PACS-ICC, Celestijnenlaan 200B, B-3001 Leuven, Belgium 
e-mail: groenSster . kuleuven . ac. be 



Received; accepted 

Abstract. Hipparcos parallax measurements of stars in the Pleiades notoriously result in a cluster distance of 
118 pc, which is approximately 10% shorter than the 'classical' result obtained from earlier main sequence (MS) 
fitting studies. In an earlier paper we developed a purely empirical MS-fitting method in an attempt to address 
this problem. This work produced conflicting results for the Pleiades between the {B — V) and {V — /) colour 
indices, indicating that the cluster's photometric metallicity is substantially lower than its generally accepted 
spectroscopic metallicity of [Fe/H] — —0.03. We were able to reconcile the discrepancy by assuming [Fe/H] — —0.4, 
the appropriate metallicity indicated from {B — V) / {V — I) colour-colour plots, and the distance moduli obtained 
from the two colour indices were in agreement with the Hipparcos result, within the la errors. 
With the release of the 2MASS All Sky Catalogue, we now apply our MS-fitting method to the Pleiades using 
the infrared colours in addition to the optical bands, in order to test the plausibility of our earlier result. 
Using the full field dwarf sample and fitting in the V/ {V — K) and K/{J ~K) colour planes, we find that assuming 
a substantially subsolar metallicity does not produce distances in agreement with the {B — V) and (V — 7) results. 
However the infrared plus (V — I) distances are in mutual agreement when adopting the spectroscopic metallicity. 
By considering only the field d warfs with My < 6, i.e. brighter than the magnitude where the Pleiades [B — V) 
colours start to be anomalous ('Stauffe r et aljl2003fl . the infrared and optical colour indices all yield consistent 
distances using the spectroscopic [Fo/H] . The concordant distances thus obtained from the V/{B — V), V/{V — I), 
V / {V — K) and K / {J — K) planes yield a mean of 133. 8±3 pc, in excellent agreement with both the pre-Hipparcos 
MS-fitting results, and the most recent determinations from other methods. 

We conclude that there are two distinct, and unrelated, issues affecting the Pleiades: 1) the Hipparcos parallax is 
in error by ~ 10%, as previously claimed; 2) the {B — V) colours of the lower MS are anomalous, and we caution 
against using the {B — V) index for MS-fitting to the Pleiades and similarly young open clusters. 

Key words, open clusters and associations: general - open clusters and associations: individual: Pleiades - 
Hertzsprung-Russell (HR) and C-M diagrams - Stars: distances 



1. Introduction 

The distance to the Pleiades h as been controversia,l ever 
since the Hipparcos mission llPerrvman fc ESAl Il997^ 
measured parallaxes of individual stars in the cluster - 
the mean distance to the cluster determined from in- 
dividual star parallaxes was found to be approximately 
10% shorter than distances previously determined from 
main sequence (MS) fitting methods. Specifically, the 
parallax measurements yield a distance of 118.34 pc, 
corre sponding to a dist ance modulus of (m — M)o = 
5.37 ijvan Leeuwenlll999l and references therein), whereas 



MS-fitting generally yields ( m — M) 



132 pc) or even longer (e.g. iM itchcU fc Johnso. 
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After the release of the Hipparcos results, several au- 
thors attempted to address this discrepancy using var- 
ious methods based on semi-empirical MS-fitting tech- 
niques, and all found distance moduli in agreement with 
the earlier MS -fitting results ('Pinsonneaul t et "all Il998l: 
ISoderblom eTal. 1998; .Stello fc Nissen 200l|) . On the 
other hand ICastellani et al. i li2002t^ . who used purely the- 
oretical isochrones in their fits, were able to retrieve the 
Hipparcos parallax distance by assuming a substantially 
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subsolar metallicity for the cluster. Concerns were raised 
that most of the distance determinations in the literature 
use methods which have some model dependence which 
may introduce unquantified systematic errors, hen ce in 
an earlier paper ijPercival. Salaris. fc Kilkennvll2003l here- 
after PSK03) we developed a purely empirical MS-fitting 
method. This method employs a sample of 54 local field 
stars, all with precise Hipparcos parallaxes, and with ho- 
mogeneous metallicity determinations on a scale consis- 
tent with that of the open clusters being studied. These 
field stars were used to construct a template MS which 
was fit to several opei i clusters, includin g the Hyades and 
the Pleiades fPSK03: IPercival fc Sal;iris..200a) . Fitting in 
both the {B — V) and {V — I)c colour planes^, this 
method precisely reproduces the Hipparcos distance to the 
Hyades, yielding (m — M)q = 3.33 ± 0.06. Applying the 
same method to the Pleiades, we found that the {B — V) 
fits yield (to - M)o = 5.76 ± 0.06 whilst {V - I) gave 
(to-M)o = 5.58 ±0.04 (corresponding to 141.9 and 130.6 
pc respectively). Not only are these results much longer 
than the Hipparcos parallax distance, they are also mutu- 
ally inconsistent. 

The location of the MS in a colour-magnitude dia- 
gram (CMD) is, of course, dependent on metallicity and 
hence MS-fitting methods rely on matching the metallic- 
ity of the template MS to that of the cluster. In PSK03 
we used the generally accepted spectroscopic metallic- 
ity for the Pleiades in our 'standard' fits, specifically 
[Fe/H] = —0.03 ± 0.06 (this val ue, plus error bar is taken 
from the catalogue of jGratton 2000 and is ba sed on the 
HRS determination of lBoesgaard fc Frie]|ll990|) . Since the 
{B — V) and {V — I) colour indices have different sensi- 
tivities to metallicity, the discrepancy found between the 
{B - V) and {V - I) MS-fitting results suggests that the 
spectroscopic metallicity is not the appropriate one to use 
when applying MS-fitting to the Pleiades. Assuming that 
both the [B — V) and {V — I) colours are normal for some 
metallicity, PSK03 used colour-colour diagrams to demon- 
strate that the photometric metallicity of the Pleiades is 
consistent with [Fe/H] = —0.4, and hence argued that 
this is the appropriate metallicity to use for MS-fits to the 
cluster. Repeating the MS-fits using this lower metallicity, 
PSK03 found that consistent distances were obtained be- 
tween the two colour planes. Furthermore, the concordant 
distances derived at this assumed metallicity are much 
shorter than those obtained from the solar abundance 
fits, and are thus consistent with the Hipparcos paral- 
lax distance (see PSK03 for full details, and summary in 
Section rOl Table^l. Hence, in PSK03 we concluded that 
our analysis did not support any mismatch between the 
MS-fitting and Hipparcos distances for the Pleiades and 
that the widely discussed discrepancy is just an artifact 
due to the cluster's {B~V) and iy—I) colours (and hence, 



^ Note that throughout this paper [V ~ I) always refers to 
{V — I)c, i.e. photometry in the Cousins system. Some clus- 
ter data used here has been converted from other photometric 
systems, as detailed in PSK03. 



photometric metallicity), which are inconsistent with the 
spectroscopic metallicity. 

However, the three most recent studies of the Pleiades, 
which use alternative distance determination methods, 
once again find long distance s in genera l agree ment with 
the earlier MS-fitting results. |Pan et all l)2004() examined 
data for Atlas (the second brightest star in the Pleiades), a 
wide binary which is resolved using optical interferometry. 
Combining orbital data with an assumed mass-luminosity 
relation, they determine a distance of 133 < D < 137 
pc, with a firm lower bound of 13 > 127 pc. The results 
are slightly model-dependent, since the masses of the two 
components must be taken from model isochrones - how- 
ever a 10% uncertainty on the mass only leads to a 3% 
uncertainty on the final derived distance, d ue to the pre- 
cision of the orbital parameters. Meanwhile. iMunari et all 
l|2004|) studied the eclipsing binary HD 23642. From ex- 
tensive new observations in Johnson B and V (to obtain 
light curves) and high resolution spectra (yielding radial 
velocities) they modelled the system and found a distance 
of 132±2 pc. These results are also slightly model depen- 
dent as the temperature of the primary star must be de- 
termined independently - this is done by comparing pho- 
tometry in many different system s (from the literature) to 
synthetic spectra. Most recently, ISoderblom et alJ (jifl^) 
have measured trigonometric parallaxes for three G and 
K dwarfs in the Pleiades, using the Fine Guidance Sensors 
on HST. Their net parallax result of 7.43±0.17 milliarcsec, 
corresponding to a d istance of 135 ib 3 pc, is in excellent 
ag;reement with both IPan et all (jifl^) and iMunari et alJ 
(|2iQM). 

Hence the problem of the Pleiades distance discrep- 
ancy still exists - whilst empirical MS-fitting using the 
(B — V) and (V — I) colours now yields a distance in agree- 
ment with the Hipparcos one (by assuming the photomet- 
ric metallicity rather than the spectroscopic one), other 
methods continue to find distances in agreement with the 
earlier, pre-Hipparcos, results. 

It should be pointed out at this stage however that 
there is a known problem with the Pleiades, which may or 
may not be related to the distance discrepancies (we hope 
to clarify this point later). This is that the lower main 
sequence stars (the K dwarfs) fall well below the position 
expected for a solar metallicity system in the V/ {B — V) 
colour-magnitude diagram. Hence these stars are either 
underluminous, or have {B — V) colours which are too 
blue for t heir (spectroscopic) m etallicity. As discussed in 
detail bv 'Stauffer et alJ l)2003|) (hereafter S03), this fact 
has been known for many years, but has largely been ig- 
nored, and certainly not explained. S03 present new spec- 
troscopic observations of several Pleiades K dwarfs and, 
by comparison with similar observations of K dwarfs in 
Praesepe, they show that the unusual blue colours in the 
Pleiades arise from anomalous spectral energy distribu- 
tions (SEDs) - at least for the two Pleiades K dwarfs that 
they study. This anomaly, which S03 ascribe to rapid stel- 
lar rotation and 'spottedness', causes the Pleiades stars 
to be approximately 10% brighter in the i3-band than 
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their Praesepe counterparts, whilst the y-band flux is un- 
affected. These spectroscopic measurements indicate that 
the {B — V) colours of the Pleiades K dwarfs should be 
about 0.1 mag bluer in {B — V) than the 'standard' solar 
sequence, in agreement with the broadband observations. 

Comparison of the shape of the MS for the two clusters 
shows that in the V/ {B — V) CMD the two sequences start 
to diverge at My ~ 6.0, in the sense that the Pleiades MS 
becomes increasingly bluer (or fainter) than the Praesepe 
MS towards fainter magnitudes (see S03, their Figure 4). 
Crucially, this divergence occurs right in the middle of the 
magnitude range of the field dwarfs used in our empirical 
MS-fitting work, which have 5.4 < My < 7.0. However, 
similar comparisons in the V/ (V—I) and V/ {V—K) planes 
show that the two sequences have identical shapes in the 
full magnitude range of our field dwarfs. This fact is very 
important for our study, as it implies that only the {B — V) 
index is affected by an anomaly, and that MS-fitting using 
(V — /) and iy — K) should give reliable (and consistent) 
results using the spectroscopic [Fe/H] for the cluster. We 
note here that when applying a MS-fitting method, for 
a fixed metallicity, the bluer (or more subluminous) the 
cluster MS, the longer the derived distance will be - hence 
the PSK03 results for (B - V) and {V - 1), at [Fe/H] = 
—0.03, are consistent with the observations of S03. 

With the release of the 2MASS All Sky Catalogue, we 
are now able to test our empirical MS-fitting method using 
the infrared colours in addition to the optical bands. We 
now know (from S03) that MS-fitting using (B—V^ colours 
may give spurious results for the Pleiades, therefore it is 
important to apply the method using the iy — K) colours 
to see whether they yield a distance which is consistent 
with the (V — I) result, since both indices appear to be un- 
affected by any colour anomalies in the magnitude range 
of interest. Also, if concordant distances are obtained, we 
vitally need to know whether they are consistent with the 
Hipparcos parallax distance, or the longer distance deter- 
mined from other methods. The 2MASS data also allows 
us to utilise the (J — K) colours. Whilst we have no spe- 
cific information on the SEDs of the Pleiades K dwarfs 
in this part of the spectrum, the (J — K) index should 
be much less sensitive to differences in metallicity than 
either {V — I) or {V — K), and is also much less affected 
by extinction. Using the K/ (J — K) CMD also ensures 
that the photometry is completely homogeneous between 
cluster and field stars and thereby minimizes systematic 
errors. Using (J — A') and {V — K) colour indices in ad- 
dition to {V — I) helps to constrain the cluster metallic- 
ity (and reddening) since the sensitivities are different for 
each colour index and consistent distances must be ob- 
tained from all the colour planes used if the results are 
valid. This should enable us to determine once and for all 
whether the discrepancy between the MS-fitting distance 
and the Hipparcos parallax distance for the Pleiades is 
real. 

As in PSK03, the method can be tested on the Hyades 
since fits using the {V — K) and (J — K) must also be 
able to reproduce the Hipparcos distance modulus for the 



cluster. We can also use M67 as a comparison cluster for 
the Pleiades since it has a similar spe ctroscopic meta llic- 
ity ([Fe/H] = 0.02 ± 0.06 according to ICrattonl EoOT?) . its 
{B — V)/{V — /) colour-colour plot is consistent with so- 
lar metallicity field dwarfs (PSK03) and its distance is not 
disputed (see Sect. l3.2|l . 

The layout of the rest of the paper is as follows: Section 
2 lists the sources of data used in this work; Section 3 gives 
a brief overview of the empirical MS-fitting method, and 
presents the results of its application to the Hyades, M67 
and the Pleiades, and Section 4 contains a discussion of 
the results and some general conclusions. 



2. Data sources 

The empirical MS-fitting method used here is described in 
detail in PSK03, to which we refer the interested reader. 
In outline, the method utilizes a sample of 54 local un- 
evolved (MS) field stars, with metallicities in the range 
-0.4 < [Fe/H] < 0.3 for which PSK03 obtained new 
BV{RI)c photoelectric photometry. All the stars in the 
sample have Hipparcos parallax measurements with er- 
rors less than 12%, Hipparcos catalogue entries were also 
carefully checked to avoid the inclusion of any binaries, 
and metallicities were determined from their Stromgren 
indices. JHK data were extracted from the 2MASS all- 
sky point source catalogue^ using a star-by-star coor- 
dinate search. Combining the 2MASS magnitudes with 
the parallax data yields the absolute magnitudes, Mj, 
Mh and Mk, which are then corrected for Lutz-Kelker 
bias, as described in PSK03. We note here the 2MASS 
JHK data have already been retrieved and utilized by 
ISaraiedini et alJ (|2004) in a MS-fitting study of several 
open clusters, and appear in their Table 2. 

2MA SS data for the Hyades were retrieved for all stars 
listed in IPerrvman et a l. (1998) as single, definite clus- 
ter members, again using a coordinate search. Data for 
the Pleiades were obtained in a similar manner, individ- 
ual MS stars having been identified from the Mermilliod's 
WEBDA dat abase^ (and see PSK03) and known bina- 
ries removed ( Bouvier et al.lll99^ iRaboud fc MermilliodI 
1998). For M67, the s i ngle s tar sequence was taken from 
Table 5 of ISandauistI l)2004|) . individual stars and their 
coordinates were then identified bv c ross-correlating with 
the data of lMontgomerv et a,l.l ill 9931) . 

Since the field star and cluster data from 2MASS are 
all in the same photometric system, for which the K filter 
used is actually Kq, [K short), no photometry conversions 
were necessary. 



3. MS-fitting using 2MASS colours 

Before MS-fitting can be performed, the [Fe/H] depen- 
dence of the colour index used must be determined so 



^ hosted by http://irsa.ipac.caltech.edu/ 

^ http:/ /obswww.unige.ch/webda/webda.html 
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that the field star template can be matched to the metal- 
licity of the cluster - here we use {V — K) and (J — K). 
This procedure is fully described in PSK03, where it was 
applied to {B — V) and {V — I), and it is also utilized by 
ISaraiedini et al.1 l)2004j) to determine the metallicity de- 
pendence of {y — K). Firstly, the slope of the relevant 
portion of the MS is estimated using the Hyades MS as 
a guide, then using this slope, each field star is 'shifted' 
along this vector to an absolute magnitude of My — 6.0. 
The colour of each star at My = 6.0 (i.e. {V — K)mv=6.o) 
is then plotted against [Fe/H] to determine the metallicity 
dependence of the colour index (see top panel of Figure^]) . 
Using this procedure we find A(T/ - JCVAfFe / H] = 0.190, 
in complete agreement with lSaraiedini et al.l l|2004|) . who 
find A(V - K)/A[Fe/R] = 0.185. 

When working with the (J — K) colour index we will 
be using the K/{J — K) CMD and hence the procedure 
is modified to find the (J — K) colour at Mk = 4.0 (this 
roughly corresponds to My = 6.0 for the field star sample 
in that it falls near the middle of the range of magni- 
tudes for the sample) . The precise metallicity dependence 
of (J — K) is harder to determine since the colour range 
of the field star sample is very small and the resultant 
(J — K)MK=i-0 vs [Fe/H] plot is dominated by intrinsic 
scatter. In fact, the total colour range of the field dwarf 
sample, at their observed colours, is only 0.25 in {J — K) 
compared with 0.91 in {V — K). 

A formal fit to the field star data yields a metal- 
hcity dependence of A(J - i4:)/A[Fe/H] = 0.078 al- 
though most of this dependence comes from the stars with 
[Fe/H] < 0.0 (see bottom panel of Figure P). Fitting the 
subsolar and supersolar metallicity ranges separately we 
find A(J - is:)/A[Fe/H] = 0.086 for [Fe/H] < 0.0, whilst 
A(J-ii:)/A[Fe/H] is virtually negligible for [Fe/H] > 0.0. 
We note here that this difference in metallicity depen- 
dence for the two metallicity regimes is in qualitative 
agreem ent w ith the theoretical i sochro nes of lGirardi et al.l 
l)2000|) '^ and IPietrinferni et alJ l)2004|) ^. However, the er- 
rors on these fits are of the same order of magnitude as 
the slopes, due to the relatively large scatter in the colours 
and, in fact, tests using the Hyades show that the metal- 
licity dependence in (J — K) is actually consistent with 
zero (see section IXT|) . 

3.1. The Hyades 

The Hyades fiducials were determined by fitting a polyno- 
mial (cubic) to all the available single star data, reddening 
was assumed to be zero and the metallicity was taken to 
be [ Fe/H] = 0.13 ± 0.06 (cluster metallicities are taken 
from lGrattoll2000l as explained in PSK03). Shifting the 
field star sample to [Fe/H] = 0.13 and fitting to the cluster 
fiducial in the V/ {V — K) plane gives a best-fit distance 
modulus of (m — M)q — 3.33 ± 0.04, where the Icr error 
accounts for photometry errors for the field stars, errors 



2.2 - 
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Fig. 1. Metallicity dependence in {V — K) and (J — K). 
Top panel - best fit relationship for the {V — K) data, 
A(y - /f)/A[Fe/H] = 0.19; bottom panel - best fit 
single slope to the (J - K) data, A(J - i4:)/A[Fe/H] = 
0.078. Note that the metallicity dependence is negligible 
for [Fe/H] > 0.0 

on magnitudes due to their parallax errors, and error due 
to the cluster metallicity, all added in quadrature. 

In the K/{J — K) plane, applying the maximum possi- 
ble metalhcity dependence (A( J - ii:)/A[Fe/H] = 0.078) 
to the field stars yields a best-fit distance modulus of 
(to — A/)o = 3.40 ± 0.04 which is marginally inconsis- 
tent with the Hipparcos parallax result of (m — M)q — 
3.33±0.01. Since the mean metallicity of the field star sam- 
ple is [Fe/H] « -0.07 and the Hyades is at [Fe/H] ^ +0.13 
(i.e. near the top end of the range of metallicities for the 
field stars), a relatively small error in the estimation of the 
metallicity dependence can become significant when the 
field stars are shifted. Hence it appears that the estimate of 
A( J-ii:)/A[Fe/H] = 0.078 is slightly too high. Repeating 
the MS-fit assuming no metallicity dependence in {J ~ K) 
(i.e. A(J - ii:)/A[Fe/H] = 0.0) results in a best-fit dis- 
tance modulus of (to — M)o — 3.35 ± 0.04 (random errors 
only), in complete agreement with the Hipparcos result. 
The uncertainty on the level of metallicity dependence in 
the (J — K) index induces a systematic error of ~ 0.06 
mag on the Hyades distance modulus since the colours of 
the stars are being 'shifted' on average 0.2 dex in metal- 
licity, with an uncertainty of 0.078 in A(J - X)/A[Fe/H], 
and the slope of the MS in the K/{J - K) CMD is ap- 
proximately 4. 

It is important to realise however that the assumed 
level of metallicity dependence makes negligible differ- 
ence to the MS-fitting results for the Pleiades. This is 
because the metallicity of the Pleiades is very close to the 
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Fig. 2. Fits to the Hyades MS. All 3 plots show the Hyades 
single star data (open triangles) and polynomial fit to the 
MS line, both shifted in magnitude by the best-fit distance 
modulus, and the field star data at their absolute magni- 
tudes (filled circles). Top panel - V/iy-K) CMD, using 
best-fit distance modulus of (m — M)o = 3.33; middle 
and bottom panels - K{J ~ K) CMDs, showing shifted 
and unshifted field stars using best-fit distance moduli of 
(m — M)o = 3.40 and (to — M)o = 3.35 respectively. 



mean of the field star sample (we assume [Fe/H]pioiados — 
-0.03 ± 0.06 in our standard fits). When building the MS 
template, individual field stars are being shifted both ways 
in the CMD, i.e. from lower to higher and from higher 
to lower metallicity, to match the metallicity of the clus- 
ter. Since the mean metallicity of the field stars is well 
matched to the cluster metallicity, most of the effects of 
the uncertainty on the metallicity dependence cancel out 
in this process, (see section . 



3.2. M67 

Fiducials for M67 were determined in the same way as 
for the Hyades, by fitting a polynomial to all the avail- 
able single star data. The cluste r metallicitv wa s taken 
to be [Fe/H] = +0.02 ± 0.06 ((Crattonl \200(t) and a 
reddening val u e of E(B — V) = 0.04 was taken from 
iTwarog et al Relative ex tinctions and r edden - 

ings were calculated according to ICardelli et al.l l)l989j) , 
so that Ak = O.lUAv, E(y - K) = 2.75E{B - V) 
and E{J - K) = 0.52E{B - V). After correcting the 
cluster fiducial for extinction and reddening, and shift- 
ing the field stars to the cluster metallicity, the best-fit 
distance modulus in the V/{V — K) plane was found to 
be (to - M)o = 9.65 ± 0.06. We note that this is in 
compl ete agreement with the results of ISaraiedini et all 
\200i\ who find (to - M)v = 9.74 ± 0.06 (correspond- 
ing to (to - M)o = 9.62, with E{B -V) = 0.04) from a 
{ V -K) MS-fit of our field dwarf sample to the M67 data 
of lMontgomerv et all (|l993). In the K/{J - K) plane the 
best-fit is (to— M)o = 9.63±0.06 if the stars are shifted us- 
ing A(J-ii:)/A[Fe/H] = 0.078, or (to-M)o = 9.61±0.06 
if no metallicity dependence is assumed. The quoted la 
errors are as for the Hyades, but now also include the ef- 
fect of an assumed uncertainty of 0.02 mag in E{B — V). 
These results are in complete agr eement with the results o f 
MS-fitting using optical CMDs -"Percival & Salaris {20Q^ 
find (to- Af)o = 9.60 ±0.09 whilst .Sandauist. (,2004) finds 
9.60 ±0.03. 

3.3. The Pleiades 

The Pleiades fiducials were determined in the same way 
as for the Hyades, by fitting a polynomial to all the single 
star data. As in PSK03, reddenings for individual stars 
were taken from iBregeJ l(l986l) (for those with no listing 
an average of E{B — V) = 0.04 was used) and Grafton's 
metallicity of [Fe/H] — —0.03 ± 0.06 was assumed in the 
standard fits. Applying reddening and extinction correc- 
tions as before and shifting the field stars to the cluster 
metallicity, the best-fit distance modulus in V/{V — K) is 
(to — A/)o = 5.67 ± 0.06, where the la error includes an 
uncertainty of 0.02 mag in E{B — V), as for M67. In the 
K/{J-K) plane, the best-fit yields (to-M)o = 5.61±0.05 
- it is important to note that in this plane, the results are 
the same, to within 0.01 mag, whether or not a metallicity 
dependence is assumed in (J — K). 

In response to the findings of S03, we returned to the 
{B — V) data and tested the effect of imposing a cut at 
Mv = 6.0, the point at which the V/{B - V) MS ap- 
pears to diverge from the 'standard' sequence. Using only 
field stars with My < 6.0 in the fit, the best-fit distance 
modulus is (to — M)q = 5.67 ± 0.06, a reduction of 0.09 
mag compared to the result from the full sample (PSK03). 
Significantly, imposing the same cut in the V/{V — I), 
V I iy — K) and K/ (J — K) planes produces the same re- 
sults as those obtained from the full sample. 
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Fig. 3. Fits to dereddened and extinction corrected M67 
MS - symbols as for Figure 12 Top panel - Vl{V - K) 
CMD, using best-fit distance modulus of (to — A/)o = 9.65; 
middle and bottom panels - K{J~K) CMD, showing 
shifted and unshifted field stars, using best-fit distance 
moduli of (to — M)o = 9.63 and 9.61, respectively. 



We recall here that PSK03 found that by assuming 
[Fe/H]p,_,,, = -0.4, the V/{B - V) and V/{V - /) 
CMDs yielded distance moduli in agreement with the 
Hipparcos parallax result (specifically, (to — M)o = 5.46 
and 5.39 respectively, using the full field dwarf sample). 
In order to test the plausibility of this cluster abundance, 
we performed the V/{V — K) and K/{J — K) fits again, 
this time using [Fe/H] = —0.4, keeping the same redden- 
ing as for the standard fits. The V/ {V — K) fit yields 
(to - M)o = 5.53 ± 0.06 whilst the K/{J - K) fits yield 
(to-M)o = 5.49±0.05 for A( J-if)/A[Fe/H] = 0.078 and 
(to — M)o = 5.61 ± 0.05 when no metallicity dependence 
is assumed (i.e. the same result as for [Fe/Hjpjg-^^g^ = 
—0.03). Imposing a cut at My = 6.0, as before, has exactly 
the same effect on the results as for the [Fe/H] = —0.03 
fits, i.e. the V/{B — V) distance modulus is reduced by 



Fig. 4. Fits to dereddened and extinction corrected 
Pleiades MS - symbols as for Figure |21 Top panel - 
V/ {V — K) CMD, using best-fit distance modulus of (to — 
M)o — 5.67; middle and bottom panels - K{J — K) 
CMD, shifted and unshifted field stars respectively, both 
using best-fit distance modulus of (to — M)o = 5.61. 

0.09 mag, to (to — M)o = 5.37, whilst the results from the 
other colour indices remain unchanged. 

All these results are summarised in Tabled 

4. Discussion and conclusions 

4.1. Summary of results 

If we accept that the K dwarfs in the Pleiades really are 
too blue in [B — V) for their spectroscopic metallicity, as 
indicated by S03, then we should not expect the {B — V) 
and {V — /) MS-fitting results to be the same if the full 
MS is used (i.e going fainter than My — 6.0). In fact, we 
should treat with caution any models which claim to find 
concordant distances in these two planes, since the colours 
of the lower MS (the K dwarfs) are obviously anomalous 
for any value of [Fe/H], whether solar or otherwise. 
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(m - M)(, 


(m - M)o 




CMD 


at [Fc/H] = -0.03 


at [Fc/H] = -0.4 


VI{B- 


V) (full sample) 


5.76 ±0.06 


5.46 ± 0.06 


VI{B- 


V) {Mv < 6.0 only) 


5.67 ±0.06 


5.37 ±0.06 


V/(V~ 


/) 


5.58 ± 0.04 


5.39 ± 0.04 


V/{V- 


K) 


5.67 ±0.06 


5.53 ±0.06 


K/{J - 


K) (with [Fe/H] dependence) 


5.61 ±0.05 


5.49 ±0.05 


K/{J- 


K) (no [Fe/H] dependence) 


5.61 ±0.05 


5.61 ±0.05 




Hipparcos result 








(van Lccuwcn 1999) 


5.37 ±0.07 



Table 1. Summary of MS- fitting results for the Pleiades. 



Excluding the [B — V) results using the whole field 
dwarf sample, the distance moduli obtained from the 
'standard' fits in all the other colour planes (i.e. V/ {V — I), 
V/{V — K) and K/{J — K)) are in agreement with each 
other within their Icr errors - and all are discrepant by 
at least 3cr with the Hipparcos result. The average of the 
{V-I), {V-K) and {J-K) fits is (to-A/)o = 5.62±0.05, 
corresponding to 133. 1±3 pc. Imposing a cut at My — 6.0 
and including the {B — V) results, the average of the fits 
in all four colour planes is (m — Af )o — 5.63 ± 0.04, or 
133.8 pc. 

Assuming [Fe/H]p;^^^^g^ = -0.4 (as we did in PSK03), 
whilst bringing the {B — V) and {V — /) results into agree- 
ment with the Hipparcos parallax distance, does not yield 
consistent results in {V-K) and {J-K). The V/{V-K) 
distance modulus is still discrepant with the Hipparcos one 
at the 2(T level, whilst the (J — K) result is slightly more 
ambiguous because of the uncertainty on the metallicity 
dependence for this colour index. However, even using the 
maximum dependence, the result of (m— M)o = 5.49±0.05 
is still discrepant by more than Icr from the Hipparcos dis- 
tance modulus. 

In summary, each of the four CMDs used here has a 
different slope on the MS and a different sensitivity to 
the effects of reddening and metallicity. We can find no 
combination of [Fe/H] and E{B — V) which produces dis- 
tance moduli which are consistent across all four colour 
planes simultaneously when the full magnitude range of 
our field dwarf sample is used. However, imposing a cut 
at Mv ~ 6.0 brings all the distance moduli into agree- 
ment when assuming the generally accepted metallicity of 
[Fe/H] = —0.03, and average reddening of E{B — V) — 
0.04. Furthermore, the average of the best-fit distance 
moduli yields a Pleiades distance of 133.8 ± 3 pc, in 
complete agreement with the results obtained from bina- 
ries in the cluster l|Pan et all 120041: iMunari et al]l2004() 
and the most recent parallax determination from HST 



4.2. Some caveats 

When applying the MS-fitting method to any cluster to 
derive its distance, using the local field stars as a template, 
we are assuming that we are comparing like with like - 



practically, this is the only sensible assumption we can 
make. The Pleiades is known to be anomalous in several 
respects - we know it is very young (~ 100 Myr) and has 
fast rotating stars, which may be affecting the {B — V) 
colours, whilst the field stars have unknown ages but are 
likely to be much older (typically a few Gyr). 

Spectroscopic metallicity determinations rely on some 
assumptions, one of which is a temperature scale. The 
spectroscopic determination used here (and generally re- 
garded as the most reliable one) is fromlBocsgaard fc Friel 
( k l990 ) , who used the temperature scale of Bochm-Vi tensa 
1)19811) . This temperature scale is based on {B — V) 
colours - i.e. from an observed colour, a temperature is 
inferred. Boesgaard & Friel determined temperatures for 
the Pleiades stars using {B — V), Stromgren h ~ y and 
H/3 photometry, all essentially in the sa me (blue) portion 
of the optical spectrum, as detailed in iBoesgaard et alJ 
l(l988l) . They themselves noted that the {R — I) data was 
also examined for use in determining the temperatures but 
go on to say that it was "rejected on the basis that all cal- 
ibrations gave consistently lower temperatures than those 
obtained using the three other indices" . We remark that 
this is consistent with the observation that the late-type 
MS stars in the Pleiades have {B — V) colours which are 
anomalously blue when compared to indices in other parts 
of the spectrum, in that a bluer colour indicates a higher 
temperature. 

The problem of determining which colour index yields 
the most appropriate temperature is not just confined to 
the Pleiades. I n a study of photo spheric abundances in 
active binaries, iMorel et al.l l|2004l) state that there is a 
tendency for the {V — R) and {V — I) indices to yield 
systematically lower temp eratures than (i? — V). As was 
the case for the Pleiades, iMorel et a,].l regard the 

{V — R) and {V — I) colour temperatures as being "spu- 
riously low". On the other hand, in a study of the effect 
of chromospher ic act ivity on the mean colour of late-type 
stars, .Amadol l|2003() find evidence for a blue excess in 
both the {U — B) and {B — V) indices for active stars 
when compared with quiescent ones of the same spectral 
type. In fact,,Amado ( 2003) cautions against the use of the 
{B — V) index when determining fundamental parameters 
(e.g. temperature) for late-type active (i.e. young) stars 
and suggests that near-infrared colours should be better 
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temperature indicators for these stars. It is beyond the 
scope of this work to predict what effect this may have on 
the derived spectroscopic metaUicity for the Pleiades and 
other young clusters, but it should be borne in mind when 
assuming a metaUicity for MS-fitting. 

Reddening estimates are also generally determined 
from a star's colours, assuming that they are normal for 
their particular spectral type. Most of the iBreeeJ l)l986tl 
reddenings for the Pleiades stars are derived from broad- 
band {B — V) or Stromgren b — y indices, both of which 
are in the potentially anomalous part of the spectrum. If 
the B flux is too strong for the late-type (lower MS) stars, 
and hence the colours too blue at a fixed spectral type, 
this would seem to indicate that the derived reddenings 
would be under-estimated. This would only serve to make 
the Pleiades problem worse in that assuming a higher red- 
dening would yield an even longer distance. However the 
referee has pointed out that the Pleiades reddening (and 
metaUicity) measurements are in fact heavily weighted to- 
wards stars of early spectral types (mostly B, A and F), 
and so the canonical reddening and metaUicity values may 
not be greatly affected by the anomaly displayed by the 
K dwarfs. 

4.3. Main conclusion 

We conclude that there are two distinct, and unrelated, 
issues pertaining to the determination of the Pleiades dis- 
tance: 

— We believe that the Hipparcos parallax distance to the 
Pleiades really is too short by approximately 10%, as 
initially suspected. 

- The {B - V) colours of the lower MS (the K dwarfs) 
are anomalously blue and we caution against their use 
for MS-fitting to the Pleiades and other young open 
clusters. 

More detailed spectroscopic study is urgently needed 
to determine the true [Fe/H] for the Pleiades - bearing in 
mind that spectroscopic metaUicity determinations have 
to assume some temperature scale, any anomalies must 
be accounted for (or at least qualitatively assessed) in this 
process. 
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